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Abstract 



. The goal of this paper is to derive a formula for the finite dimensional joint characteristic 

p I ' function (the Fourier transform of the finite dimensional distribution) of the coupled process 

{(Wt,Lf) : t £ [0, oo)}, where {Wt : t € [0, oo)} is a d-dimensional Brownian motion and 
{Lf : t £ [0, oo)} is the generalized d-dimensional Levy's stochastic area process associated to 
a d x d matrix A. Here A need not be skew-symmetric, and in our computation we allow A to 
vary. The problem finally reduces to the solution of a recursive system of symmetric matrix 
Riccati equations and a system of independent first order linear matrix ODEs. As an example, 
the two dimensional Levy's stochastic area process is studied in detail. 

^" : 1 Introduction 

<N 

The Levy's stochastic area process {Lt : t G [0, +oo)} associated with a two dimensional Brownian 
motion {(W t {1) , W t {2) ) : t £ [0,+oo)} was introduced by P. Levy in [10] defined as 

(N 



L t := \ [\wPdWP - W^dWP), t > 0. 
2 Jo 



The geometric meaning of the process Lt is the algebraic area enclosed by the two dimensional 
Brownian path up to time t and the segment of the origin O and the point (W t , W t )• By 
using stochastic Fourier expansion of Brownian motions, Levy derived the conditional characteristic 
function of Lt with respect to Wt as 



E[exp(i\L t )\W t = x]= . > exp 

2 sinn(iA/2j 



\x\ 2 f tX tX 

— 1 coth — 

2t \ 2 2 



It follows that the characteristic function of Lt is given by 

and the joint characteristic function of the coupled process {(Wt,Lt) : t ^ 0} can be computed 
explicitly. 
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To compute the finite dimensional joint characteristic function of the coupled process is harder 
than the one-dimensional marginal one, since it involves complicated correlations among different 
time spots; although the process is Markovian, it is still nontrivial to derive the transition density 
explicitly. 

In T. Hida's paper [5], he considered the problem from a different and a more general perspec- 
tive, under the framework of Wiener-Ito's chaos decomposition. He studied the law of quadratic 
functionals of Brownian motion on the canonical Wiener space W, which are elements in the sec- 
ond order chaos component. The classical Levy's stochastic area process is a special case in his 
general setting. The main idea of his approach is to use the celebrated Wiener-Ito's chaos de- 
composition theorem to construct a one-to-one correspondance between such quadratic functionals 
and symmetric Hilbert-Schmidt operators on the Cameron-Martin subspace H of W . By a formal 
computation, he found out that for any quadratic functional F on W, if B is the corresponding 
symmetric Hilbert-Schmidt operator on H, then the characteristic function (Fourier transform) of 
F is given by 

oo 

/ e <AF ( w V(dw) = TT( 1 " 2i\X n y 1/2 e- lXXn , 

J W n=l 

where {A n : n ^ 1} is the family of eigenvalues of B with multiplicity. He also pointed out that the 
right hand side of the above identity is equal to 

(det 2 (/-2iA J B))- 1 / 2 , 

where det 2 denotes the regularized Fredholm determinant. By studying the special example of 
Levy's stochastic area process and computing the eigenvalues of the corresponding operator, he 
recovered Levy's formula. 

After Hida's important work, different methods of computing the regularized Fredholm deter- 
minant in order to study the laws of a wider class of quadratic Wiener functionals were developed 
by several Japanese mathematicians: N. Ikeda, S. Manabe, S. Kusuoka, K. Hara, etc.. The funda- 
mental ideas in their works can be summerized as two approaches: to compute the spectrum of the 
corresponding Hilbert-Schmidt operator directly, or to reduce the infinite dimensional case to the 
finite dimensional one by posing additional assumptions on the operator. The first approach was 
developed in Ikeda and Manabe's paper [7], in which they computed the spectrum of a variety of 
quadratic Wiener functionals in order to study the asymptotic behavior of stochastic oscillatory in- 
tegrals on the Wiener space. The second approach was developed from different aspects based on a 
fundamental decomposition assumption proposed by Ikeda, Kusuoka and Manabe. They restricted 
themselves to the study of operators which could be divided into an operator of finite dimensional 
range and an operator of Volterra's type, in order to reduce the case to a finite dimensional one 
in a certain sense. In Ikeda, Kusuoka and Manabe's paper [6], based on such decomposition, the 
computation of the regularized determinant is reduced to the computation of a finite dimensional 
determinant, in which algebraic methods and differential equations methods can be applied. Later 
in [8] they developed a general method to compute the law of the corresponding type of quadratic 
Wiener functionls by using the Cameron-Martin transformation along critical paths. In the mean- 
while, they put forward another idea for the computation by relating the problem to the Van 
Vleck-Pauli formula in quantum mechanics and derived a formula for the law of a certain class of 
quadratic Wiener functionals by using physical approach. In the paper [3] of Hara and Ikeda, they 
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further developed the ideas in [8j by relating the problem to the study of dynamics on Grassman- 
nians. Due to the fact that the determinant is expressed by the solutions of the Jacobi equation 
with Van Vleck-Pauli's formula, they showed that the determinant could be computed in terms of 
Pliichker coordinates of Grassmannian manifolds. 

From the works on computing the regularized determinant, one can see that the computation is 
very involved, even reduced to the finite dimensional case. In the series of works mentioned based 
on Hida's formula, they considered the marginal law of the two dimensional Levy's stochastic area 
process and recovered Levy's formula. 

It should be pointed out that after Levy's original work, several mathematicians proposed 
different methods to compute the marginal characteristic function of the two dimensional Levy's 
stochastic area process directly. In particular, B. Gaveau [2] studied the marginal distribution of 
the process (inversion of the characteristic function). In the present paper, we are going to study 
the finite dimensional joint characteristic funtion of the d-dimensional generalized Levy's stochastic 
area processes from a different angle of view by computing the multi-dimensional Fourier transform 
directly. Our work is based on the idea of K. Helmes and A. Schwane in their paper [3] for the 
computation of the marginal characteristic function of the d-dimensional Brownian motion together 
with the general d-dimensional Levy's stochastic area processes. 

Let {Wt : t G [0, oo)} be a d-dimensional Brownian motion, and let so(d) be the space of d x d 
skew-symmetric matrices, where d ^ 2. Fix A G so(d), the process 



is called the general d-dimensional Levy's stochastic area process associated to A. In Helmes 
and Schwane's paper [4], they derived a formula for the marginal characteristic funcion of the 
coupled process {(Wt,Lf-) : t G [0, oo)}, based on the Cameron-Martin-Girsanov's transformation 
theorem and Ito's formula. In our paper, we are going to further exploit this method to establish 
the finite dimensional joint characteristic function (the Fourier transform of the finite dimensional 
distribution) of the coupled process in a more general setting. We will see that it reduces to the 
solution of a recursive system of symmetric matrix Riccati equations and a system of independent 
first order linear matrix ODEs. 

Since the computation is quite lengthy, we first present the main result in our paper and the 
main idea for the proof. All details and technical steps are left to section 2. As an example, we 
study the two dimensional Levy's stochastic area process in detail. 

Assume that {Wt,J~t '■ t 0} is a d-dimensional Brownian motion on some probability space 
(f2,.F, P). For any A being a real d x d matrix (not necessarily skew-symmetric), we define the 
generalized Levy's stochastic area process {Lf : t ^ 0} associated to A as 



J o 

Fix = to < t\ < t2 < ■ ■ ■ < t n < oo, n ^ 1. Assume that Ax,'" , A n are real d x d matrices. For 
7li " ' j In ^> A.l> • • • , A n G K, we are going to compute the following functional: 





/(7 
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where i = \f— T. Notice that the matrix can vary on different time spots. 

Throughout this paper, < h < ■ ■ ■ < t n < oo and A u ■ ■ ■ , A n G M d (R) wiil be fixed. The 
notation * wiil denote the transpose operator. We should point out that even for complex matrices, 
* is simply the transpose without taking conjugate. For z = (z 1 , • • • , z 3 ) G C™, we use (z) 2 to denote 
Y^j=i( z *) 2 which is the holomorphic extension of the real case. TV will be denoted as the trace 
operator. 

Our main result is the following. 

Theorem 1.1 The functional f("fi, • • • , 7 n ; Ai, ■ ■ ■ , A„) is determined by 



fdi,--- ,7 n ;Ai,--- ,A n ) 



f[expd f tj Tr(K iAj ,..., iAn (s))ds 

3=1 ° 

~\ T \,' -- ./A, '• s '''"/.\ . •• ./A, (tj ''/'./ ~d*\. 



Here {Ki Aj> ... ,i An (t) : t G [0, = 1, • • • ,n} is defined as the solution of the recursive system of n 
symmetric matrix Riccati equations (starting from j = n to j = I): 



d 

■^K iAjt ..., lAn (t) = C iAj ,...,iA n {t) ~ K iAj ,...,iA n (t){ ^2 K iK,- ,iA n {t) 

r=j+l 

n n 

+ Y^(iKA r ))-{Y, K iAr ,.., iAn (t) 

r=j r=j+l 
n 



+ J2(^rA* r ))K iAj ,..., iAn (t) - Kf Ajj ... jiAn (t), 



r=j 



t g [o,tj], 

K iAj ,...,i An (tj) = 0, 



where 



C iAj ,.., iAn (t) := A?AjAj — iAj[( £ K iAr ,..., iAn (t) + £ {iKKMj 

r=j+l r=j+l 
n n 

+A*( £ K tAr> ... yiAn (t) + ( iA rA))], t G [0,tj]. 

r=j+l r=j+l 

Moreover, {Hi Ajt ... t i An (t) : t G [0,tj],j = 1, ■ ■ ■ ,n} is i/ie solution of the system of n independent 
linear matrix ODEs: 

d n n 

J t H iAj ,.., iAn {t) = {J2 K iA r ,-..,iA n (t) + Y^{ iA rA r ))H iAj ... :iAn (t), 
r=j r=j 

t€ [0,tj], 
H iAj ,.., iAn (0) = Id, 

3 = !,■■■ ,n, 
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and {fij : j = 1, • • • , n} is defined recursively as 

j = n-l,-- - ,1. 

The subscripts of -RT and are used to indicate the dependence on A. We do not wish to 
address the existence and uniqueness of the preceeding maxtrix differential system of Riccati type. 
Under our context, there is a unique solution to the system appearing in Theorem ll.il for details, 
the reader may consult [9], [13], [H] and etc. 

Now we are going to briefly express the main idea of the proof in an informal way. 

The fundamental tool of the proof is Girsanov's transformation. However, it fits for real-valued 
processes only. Hence we focus on the real case first, that is, consider the functinal 

n n 

dill,-- - ,-y n ]Xx,--- ,A n ) := E[exp{^i( 7fc , W tk ) + ^A fc L t ^ fc }] 

k=l k=l 

first for small Ai,--- , A n G R, and then try to complexify the result by standard arguments in 
complex analysis. 

Our idea of computing the functional g is to eliminate the stochastic integrals one by one, 
starting from the largest time interval [0, t n ], by using the Girsanov's transformation theorem, and 
try to track the original Brownian motion along each time of transformation in order to handle the 
W . When applying change of measure, a stochastic integral (with respect to a proper Brownian 
motion) is transformed to one-half of its quadratic variation, which is a Lebesgue integral almost 
surely. To handle this term, we introduce a symmetric matrix Riccati equation to split it into three 
terms, which will be transformed to a deterministic one by change of measure again and applying 
Ito's formula for a suitable process. Recursively, we will be able to cancel out all of the n stochastic 
integrals in g, and obtain a recursive system of symmetric matrix Riccati equations. To handle the 
W term, we need to track the change of W along every time of transformation. From Girsanov's 
theorem, we will see that the diffusion form of W under each transformation is actually invariant 
(it is always Gaussian by solving a linear SDE), which will enable us to do the computation easily. 
This procedure will lead us to a system of independent first order linear matrix ODEs. 

Our interest in looking for an (as far as possible) explicit formula for the joint law of Brownian 
motion together with its Levy area is motivated by the recent progress in the understanding of the 
solutions of Ito's stochastic differential equations revealed recently in T. Lyons' work |12| . in which 
it has been demonstrated that a large class of Wiener functionals (including Ito solutions to SDEs) 
are continuous with respect to Brownian sample paths and their Levy areas equipped with the law 
of Brownian motion together with the Levy area process. For details about these developments, 
see T. Lyons [12\, T. Lyons and Z. Qian [11], P.K. Friz and N.B. Victoir [J. 

2 Proof of the Main Result 

Now we are going to work out the idea in section 1. The procedure is to handle the real case first 
and then use complexification. 

The following technical lemma is important (see also [1]), which gives us finiteness in the real 
case. 
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Lemma 2.1 There exists a number c > 0, such that 

n 

sup{E[exp{^ XkLf k k }} : Xj G (-c,c),j = 1,2,- •• ,n} < oo. 
fe=i 

Proof. It suffices to consider the case when n = 1. By Cauchy-Schwarz's inequality, we have 

E[exp{AL^}] = E[exp{A f {AW S , dW s }}} 

Jo 

= E[exp{A f (AW s ,dW s ) - A 2 f \AW s \ 2 ds} + A 2 [ \AW s \ 2 ds}] 
Jo Jo Jo 

sC E 1 / 2 [exp{2A / (AW s ,dW s )-^- [ \AW s \ 2 ds}] ■ 
Jo 1 Jo 



E 1 / 2 [exp{2A 2 / \AW s \ 2 ds}}. 
Jo 

Notice that the first term inside the expectation E on the right hand side is a local martingle, the 
finiteness of E[exp{AL^}] will follow immediately once we show that 

E[exp{A 2 f \AW s \ 2 ds}} < oo 
Jo 

when A is small enough. By Jensen's inequality, it remains to show that 

sup E[exp{A 2 |^W s | 2 }] < oo 
se[o,t] 

when A is small enough, which is obvious by simple calculation based on Gaussian random variables. 
■ 

Now consider small Ai, • • • , A n G R as in lemma 1, 7 l5 ■ ■ ■ ,7 n el, and the functional 

n n 

gilir-- ,7„;Ai,-" ,A n ) =E[cxp{J2i(lk,W tk } + J2^kLf k k }} 

k=l k=l 

defined in section one. The following proposition gives the formula of g. 
Proposition 2.2 The function 5(71, • • • , 7 n ; Ai, ■ ■ ■ , A n ) is determined by 
g(ii,--- ,7«;Ai,--- ,A n ) 

= fW| f 3 Tr{K X] ,.. M {s))ds-\f J // A , l .. A ,i>)//A,...A,f,i/'/ J ' / -}- 

j=l Jo Jt i-i 

Here {K x iXn (t) : t G [0,tj],j = 1,2, ■■■ ,n} is defined recursively (starting from j = n) by the 
symmetric matrix Riccati equation 

, n n 

-K Xj ,... M {t) = C Xj ,.., Xn (t)-K Xj ,... An (t)(J2 K Xr ,.., Xn (t) + Y,(KA r )) 

r=j+l r=j 

n n 

r=j+l r=j 
t€[0,tj], 
K Xj ,..., Xn (t j ) = 0, 



where 

n n 

Ca,,.,a„W : = -X^Aj-X^ K Xr ,.., Xn (t)+ (*rK))Aj 

r=j+l r=j+l 
n n 

+A*( £ K Xr ,..., Xn (t) + £ (X r A r ))}, t G [0,t 3 ] 

r=j+l r=j+l 

is also defined recursively starting from j = n. Moreover, {H x . ... Xn (t) : t £ [0,tj],j = 1,2, • • • ,n} 
is the solution of the system of n independent first order linear matrix ODEs 

, n n 

-H x .,.. M (t) = i»A.,.,A..!n • ^(A ; ,1 ; ;)//. V ..... A .. t€ [0,tj], 

r=j r=j 

Fa,,..,a„(0) = Id, 

j = 1,2, ■■■ ,n, 

and {[ij : j = 1, 2, ■ ■ ■ , n} is defined recursively by 

■■= Ij + #a7+i,- ,A„fe) iJ A j+1 ,... ,A„ (tj+l)li j+ i, j = 1, 2, ■ • • ,n - 1. 

Proof. We divide our proof into two steps. 
(1) Step one. 

Consider first that 7i = ■ ■ ■ 7 n = 0, and write 

n 

h{Xi, • • • , A n ) := E[exp{^A fc L t f }]. 

fc=i 

= E[exp{VA fc /"'VfcW^dW,)}]. 

By changing the original probability measure P on the largest time interval [0, t n ], we have 

/i(Ai,--- ,A n ) = E[exp{VA fc / (A k W s , dW s ) + X n / (A.W^dW,) 

fc=i Jo Jo 

\2 />t n \2 rt n 

-fJ o \A n W a \ 2 ds + -£ J \A n W s \ 2 ds}} 

"zi /•** a 2 /•*» 

= E n [exp{VA fe / (A k W 8 ,dW a ) + -± / |A„W s | 2 d S }], 

where E n denotes the expectation under the probability measure 

dP n := exp{A n / (AnW^dW,) - -f / | A^| 2 d S }ciP. 
jo 1 Jo 
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Notice that {Wt : t G [0, t n ]} may not be a Brownian motion under the new measure P n . However, 
by the Girsanov's theorem, under P n , the process 

W t (n) := W t - X n f A n W s ds, t G [0, t n ], 
Jo 

is a Brownian motion, and the original process {Wt : t G [0, t n ]} satisfies the following SDE: 

dW t = X n A n W t dt + dW t (n) , t G [0, i n ] . 

In order to eliminate the integral over [0,i n ], let C n (t) := —X n A n A n (t G [0, t n ]) and introduce the 
following matrix Riccati equation: 

j f Kx n (t) = C n (t)-X n [K Xn (t)A n + A* n K Xn (t)]-Kl(t), t € [0,t n ], 

K Xn (tn) = 0. 

From now on, to simplify our notation, we will use K n to denote K\ n , and later we will also 
use Kj to denote K\ jt ... t \ n , j = 1, ■ ■ ■ , n — 1 to omit the indication on the dependence on A. By 
symmetricity of the equation, the unique solution {K n {t) : t G [0,t n ]} is symmetric. Hence, by 
substitution, we have 



h(Xi, • • • , A n ) — 



E„[exp{ V A fe / k (A k W s , dW s ) - \ f " W* s C n (s)W s ds}\ 
fr[ Jo 2 J 

E n [exp{^A fc y o (A k W s ,dW s )--J^ W* s —K n {s)W s ds 

~\t W ^ K n( s )( X nA n ) + X n A* n K n (s)]W s ds 
-\J^ \K n (s)W s \ 2 ds}] 

E„[exp{^A fe y o ( Ak W s ,dW s )--J^ W* s -K n {s)W s ds 

-/ (K n (s)W s ,X n A n W s )ds-- / |^ n ( s )H/,| 2 d S }]. 
jo 2 Jo 
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By changing of measure again, 

'«* „„. 1 d 



'JL^ ftk | ftn d 

h(\i,--- ,X n ) = E n [eMYl X kj o (A k W s ,dW s )--J^ W* s —K n {s)W s ds 

- [ tn (K n (s)W s ,X n A n W s )ds- f tn (K n (s)W s ,dW^}}} 
Jo Jo 

™zi r*k i /•*« w 

= E n [exp{^A fc y o (A k W s ,dW s )--J^ w;-K n (s)W s ds 



k=l 

in 

(K n {s)W s ,X n A n W s )ds 



(K n {s)W s , dW s - X n A n W s ds)}\ 

rt k 



= E„[exp{VA fe / (A k W s ,dW s )-- W;-K n (s)W s ds 
k=1 Jo 1 Jo as 

rt n 

- / (K n (s)W s ,dW s )}]. 
Jo 

where E n denotes the expectation under the probability measure 

rtn i rt n 

dP n := exp{ / (K n (s)W s , dWjp) - - / \K n (s)W s \ 2 ds}dP n . 
Jo 1 JO 

Under P n , the process 

W t {n) := W t {n) - ( K n (s)W s ds, t G [0,t n ] 
Jo 

is a Brownian motion, and the original process {Wt : t G [0,i n ]} satisfies the following SDE: 

dW t = (K n (t) + A n ^ n )W t dt + dW t {n \ t G [0, t n ]. 

It should be pointed out that under P n , the quadratic variation process of the semi- martingale 

{Wt : t G [0, t n ]} is actually the same as that of the Browinian motion {w/™^ : t G [0, t n ]}. Now let 
F(t, w) : [0, t n ] xl^Rbe defined as 

F(t,w) := w*K n (t)w, 

by applying Ito's formula to the process {F(t,Wt) : t G [0,t n ]}, we have 

rt n j ft™ 
Jo ^^^W^ + 2 / (^«( s )^' ^ + / Tr(K{s))ds = 0, 

where Tr denotes the trace operator. Therefore, we arrive at 

h(X lr -- ,A n ) = exp{i I'" Tr(K n (s))ds} -Eniexpiy^Xk f" {A k W s ,dW s )}\. 
2 Jo k=i Jo 



Now we proceed a similar argument over the second largest time interval [0, t n -\\. The difference 
here is that {Wt : t G [0,i n _i]} is not a Brownian motion under the probability measure P n . 
However, still by changing of measure, we have 

1 f tn 

h(Xi, • • • , A„)exp{-- / Tr(K n (s))ds} 
1 Jo 

=^[exp{VA fc / (A k W s ,dW s ) + \ n - 1 / (Vi^,<^)}] 

n—2 „^ fe rt n —i 

=£„[exp{ V A fc / (A fe W s , dW s ) + A n _! / (4_iW s , + \ n A n )W s )ds 

k=i Jo Jo 



Jo 1 Jo 

\2 rtn-i 

+ |A»-iW s | 2 cfe}] 

2 JO 

n—2 „i k rtn—i 

=E„_ 1 [exp{V A fc / (A k W s ,dW s ) + A n _i / (A n ^W s , {K n {s) + A^V^ds 

+ |4»-iW s | 2 cfc}]. 

2 JO 



v 2 . rtn-i 
'o 

Here E n _i is the expectation under the probability measure 



W t (n_1) := W t (n) 



/•*n-l "~T~i A 2 /"*n-l 

dP n _! := exp{A n _! / (i^rf)-^ |A t _!^| 2 d S }(iP n . 
io 2 Jo 

Under P„_i, the process 

An_i / A n ^W s ds, tG[0,t n _i] 

JO 

is a Brownian motion, and the process {Wt : i G [0,i„_i]} satisfies the following SDE: 

dWi = (tf n (t) + A n .4„ + Xn^A^Wfdt + ^ (n ~ 1} , t G [0,t n -l]. 

Let 

C„_i(t) := -^-i^-iVi - A„_i[(K„(t) + A n ^;)^l„_i 
+A* n _ 1 (K n (t) + \ n A n )}, t G [0,tn_i], 

and introduce the following matrix Riccati equation: 

j t K n -i{t) = C„_i(t) - #„_!(*)(#„(*) + A„A n + A„_i^ n _i) 

-(K„(t) + A n < + A n _i<_ 1 ) J PC n _i(t)-K 2 _ 1 (t), te [0,t n _i] 
^„_i(t n _i) = 0, 
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we have 



fc(Ai, ■ ■ ■ , A„) exp{-i jf " Tr(K n (s))ds} 

=E„_ 1 [exp{VA fe / (AfcW^.dW*)-- / W;C n _i(s)W s <fe}] 
k=i Jo 2 Jo 

=E n _ 1 [exp{^A fe y o (AfcW,,^)--^ W;-^ n _!( S )W s cZ S 

ptn-l 

- / (^ n _i(s)W s ,(K n (s) + A n ^ n + A n _i^ n _i)W s )ds 

JO 

1 ftn-l 

--J \K n ^{s)W s \ 2 ds}] 

^ n-2 tjt /-tn-i J 



tn-1 

<iir n _i(a)W 8 , + \ n A n + An-i^iJW^da 

o 

* n " 1 (if n _ 1 ( S )W s ,cZ^- 1 ))}] 



______ «____2 -t fc .. nt n -i j 

=S n _ 1 [exp{^ A fc ^ (,W S , dW a > - - ^ W;— ^.iCaJWads 

rtn-l 

- / (ViW^,,^.)}]. 

JO 

Here we have changed the probability measure from P n -\ to 

dP^i :=exp{ ^ Vn-iW^,,^ 1 ') ~ \ T 1 \K n ^{s)W s \ 2 ds}dP n ^. 
Jo 1 Jo 

By applying Ito's formula to the process {W^ K n -i(t)Wt ■ t £ [0,t n -i]} and noticing that the 
quadratic variation process of the semi-martingale {Wt : t G [0,t n _i]} is the same as that of a 
Brownian motion, we again have 

ftn-l W ftn-1 ftn-l 

J W* s —K n ^{s)W s ds + 2 J (K n ^(s)W s , dW s ) + y Tr(K n ^(s))ds = 0. 
Therefore, we arrive at 

fc(Ai, • • • , A n ) exp{-^ jf " Tr(K n (s))ds} 

= exp{- / Tr(K n _ 1 ( S ))d S }-E n _ 1 [exp{VA fc / <A fc W 8 ,<iW 8 >}]. 
IJo t^i Jo 
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Now the recursion is quite obvious from the key observation that when applying transformation of 
probability measures, the original process {Wt}(over the proper time interval) is always a diffusion 
of the same kind, namely, it satisfies an SDE of the form 

dW t = Q(t)W t dt + dB u 

where {B t } is a Brownian motion under the corresponding probability measure. To be more precise, 
after j steps, we will have a system of j matrix Riccati equations defined recursively, 

C n {t) := -X 2 n A* n A n , t€[0,t n ], 

f t K n {t) = C n (t) - X n [K n (t)A n + A* n K n (t)} - Kl(t), t e [0,t n ], 
K n {t n ) = 0; 

'C„_i(t) := -Xl^A^An^ - X n ^[{K n {t) + A n ^*)^ n _i + A* n ^{K n {t) 

+X n A n )}, t€[0,tn-i], 

< f t K n ^{t) = C n - 1 {t)-K n _ 1 (t){K n (t) + X n A n + X n - 1 A n _ 1 ) 

-(K n (t) + X n A* n + K-iK-i) K n-i(t) - Kl^t), t e [0,t n _i], 
/n-ife-i) =0; 



C n -j+i(t) 



-X 



n—j+1 



An-j+l-An-j+l - ^n-j+l[(J2r=n-j+2^r(t) 



A 
dt J 



+ Y^r=n-j+2^rA* r ))A n ^j + i + A^_ j+1 (Y^ =n _j +2 K r (t) 

+ E"=„-, +2 (V4r))], t€[0, 

T K n _ j+l {t) = C n - j+1 (t) - K n _ j+l {t){Yr r=n - j+ 2 Kr(t) + Er=n- j+ l(KA r )) 

-(E"=„- J+ 2 K r (t) + Er=n- j+ l(KA* r ))K n _ j+1 (t) - K*_ j+1 (t), 

t e [o, t n -j+i], 

/n-j+l(in-3+l) =0, 
and we will arrive at 



h(X 1 ,--- ,X n ) 



:ex P<9 E / Tr(K r (s))ds}-E n _ j+1 [exp{J2 X * / (A k W s ,dW s )}}. 



k=i 



r=n-j+l " u 

Here under the probability measure P n _j + i, {W t '■ t € [0, t n -j + i]} satisfies the SDE 



dW t = { K r(t)+ (KA r ))W t dt + dwl n - j+l \ t e [0,t n - j+1 ], 

r=n—j+l r=n—j+l 
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where {W t (n : t G [0,t n _ j+ i]} is a Brownian motion under P n _ J+ i. By carrying out a similar 
argument, that is, by changing measures and applying Ito's formula, we will obtain 



1 n f tr 
h(X lr -- ,A„)exp{-- Yl / rr(tf r («))ds} 

= exp{- / Tr(K„_ j ( S ))d S }E ri „ j [exp{ V X k / (A fc W s , dW a )}], 
2Jo k=i Jo 



where {K n _j(t) : t G [0, i n _j]} is the solution of the matrix Riccati equation 

d n n 

jKn-.it) = C n .j{t)-K n .j{t){ K r {t)+ Y^rA r )) 

r=n—j+l r=n—j 
n n 

-( R r( t )+ E (KKWn-jW-K^it), 

r=n—j+l r=n—j 
t G [0,-tn-j], 

K-n—j{tn—j) = 0) 

in which 

n n 

C n -j(t):= -X^A^A^j - X n ^[( Y R r(t)+ E i^rA*))A n -j 

r=n—j+l r=n— j+1 

n n 

+A* n _ j ( Y R r(t)+ E i^rA r ))], t € [O^-j]. 
r=n—j+l r=n—j+l 

By induction on j, the proof of the case where 7 1 = • • • = 7„ = is now complete. 
(2) Step two. 

Now we consider the case with 7 1; • • • ,7 n 6 I. In step one, the ultimate goal of applying those 
transformations of probability measures is to eliminate the stochastic integrals one by one, starting 
from the largest time interval. After each transformation, the distribution of Wt k (k = 1, • • • ,n) is 
changed. In order to work out the case involving r y 1 , • • • , 7 n , we need to track the original process 
{W t : t G [0,t n ]} after each transformation. The main difficulty comes from the fact that if we 
apply a transformation on [0, ifc], the distribution of {Wt} over [tk, t n ] is hard to compute. However, 
by using the crucial observation that the diffusion form of {W t } is invariant, w can factor out the 
term over [tk , t n ] . 

Let's formulate the idea in detail. By using the same notation as in step one, we have 

dill,-- - ,ln^l^n) 
n 

=E n [exp{Y,i("fk,W tk )} ■ A n ], 
k=i 

where 

A n := exp{i f " Tr(K n (s))ds + VaJ " (A k W s , dW s )}. 
1 Jo k=i Jo 
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Under P n , the process {Wt : t G [0,i n ]} is a diffusion of the form 

dW t = (K n (t) + A n ^ n )W t dt + dwi n \ t e [0, t n }. 
Let {-f^AnW : * ^ [0>*n]} be the solution of the following linear matrix ODE 

j t H Xn {t) = (K n (t) + X n A n )H Xn (t), t€[0,t n ], 
Hx n (0) = Id. 

Then by the explicit solution of linear SDE, we have 



W t = H Xn (t) f H-^s)dwl; n \ t e [0,-Q. 

J 



Hence 



gill,-" ,7n^l,^n) 



=E n [exp{^ i( 7fc) W tfc > + i( ln , H x Jt n ) / fT^( S )dwi n) )} • A n ] 
k=i Jo 

=E n [exp{J2i(l k ,W tk ) + i( ln ,H Xn (t n ) / H^(s)dWr) 
k=i Jo 

+ i( ln ,H Xn (t n ) P H^(s)dW^)} ■ A n }. 

Jt-n-l 

Notice that the stochastic integral f* n i H^(s)dWs n ^ is independent of the rest since the integrand 
is deterministic, we have 



:E n [e W {i(ln,H Xn (tn) P H^(s)dW^)}]- 

Jtn-l 

n-1 

E n [exp{^( 7jfc , W" tfc ) +i(H* x ; 1 (t n - 1 )H* Xn (t n ) ln ,W tn _ 1 )} ■ A n ] 



k=i 

n-2 



~-R n • E n [exp{^ i( lk , W tk ) + <Ai n _!, Wt^)} ■ A n ], 



k=i 



where 



R n := E n [exp{i( 7n)J ff A „(t n ) / " H^(s)dW^)}}, 

Jtn-l 

A*n-1 : = 7n-l + - ff C 1 ( t ™- 1 ) F A Il ( t «)7n- 

Now we can see that the random term over [i n _it n ] is factored out. 
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Similarly, by applying transformations as in step one, we further have 
g(ii,--- ,7 n ;Ai,--- ,A n ) 



n-2 



--R n • E n _ 1 [exp{ J] i( 7fc , W tfc ) + (/x n _ 1; W^)} • Vi], 



k=i 



where 



I ftn-l 1 ftn 

A n _i: = exp{-j[ Tr(K n -i(s))d8 + - J Tr{K n {s))ds 
+ VA fe / (^ fe W s ,dW s )}. 
By a similar argument, let {H\ n _ 1 ,\ n (t) : t G [0,i„„i]} be the solution of the equation 

t € [0,t n -i], 
^-i,a„(0) = /d, 

we obtain that 

ghi,--- ,7«;Ai,-- - ,A n ) 

n— 3 

— Rn ' Rfi—1 

■ E„_i[exp{^z(7 fc , W tk ) + i(fi n _ 2 , W tn _ 2 )} ■ A n _i], 
fe=i 

where 

Jtn-1 

Finally, by a simple induction argument, we arrive at 

5 ( 7l ,--- , 7 „;Ai,--- ,A„) = n(^-exp{- / Tr(tf»)<fo}). 

j'=i 

Here 

R 3 : = Z J oK V {r llj .IIy..., x jl j] f ll x l.. x is),l\V: j }. ; = l,2,-..,n, 
{i^Aj •••,A„(i) : * G [0 , t j ] } is the solution of the equation 

A n n 

-H x .... M (t) = &K j (t) + £(\ j A j ))H x ...., Xn (t), t€[Q,tj], 

r=j r=j 

H x ...., Xn (0) = Id, 
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and {fj,j : j = 1, 2, ■ ■ ■ , n} is denned recursively by 

Mj : = 7j + An (tj)H\ j+u - ,a„ J = 1,2, • • • ,n - 1. 

It remains to compute = 1, 2, ■ ■ ■ , n) explicitly, which is easy since everything here is Gaussian. 
Namely, we have 

Rj = EjieMiiH^.^xMnj, T H'^^dW^}} 

= %[exp{z ^ (Hl-^^mi^.^^dW^)}} 



exp {-i [ tj //;,.!.., v ,i>-i// A ,,A.iO!/'r'M- 



Therefore, the proof is now complete and we have 
g(ii,--- ,7™;^!,--- ,A n ) 



= IIexp{i Tr{K 3 ( S ))d S -\ f \H^„ M (s)H* Xj ,.. ^Jt^fds}. 
j= i z J o z J tj-i 

m 

From the proof of the above proposition, we can see that the computation of 

n n 
k=l k=l 

for small Ai, ■ ■ ■ , X n € R reduces to the solution of a recursive system of symmetric matrix Riccati 
equations and the solution of a system of independent first order linear matrix ODEs. If j 1 = 
■ ■ ■ = 7 n = 0, then we don't need the ODE system at all. Now we are going to complexify our case 
before. 

Lemma 2.3 Fix j 1 , ■ ■ ■ ,7 n £l. Then when c is small enough, the function 

n n 

,z n ) :=E[e W {J2i(lk,W tk ) + Y,z k L? k k }} 
k=l k=l 

is holomorphic in the domain D c := {(z±,--- ,z n ) G C™ : Re(zj) £ (— c,c),j = 1,2, •■• ,n} ofC n . 
Moreover, the function 



V>(Ai,--- ,A n ) : 



[]exp{- / Tr{Kj(a))d8 

3=1 J ° 



defined on W 1 can be extended holomorphically toC n . Such an extension is unique, and when 
restricted to D c , 
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Proof. By lemma 1, when c is small enough, (j)(z\, ■ ■ ■ , z n ) is well defined on D c . The continuity 
of 6( zi, ■ ■ ■ , Zn) follows easily from uniform integrability. Moreover, since the function 

n n 

(*!,•••,*»)'-»• exp{^i( 7fc , W tk ) + ^2z k Lf k k } 

k=l k=l 

is holomorphic on C n for every uj G f2, by Fubini's theorem and Morera's theorem, <j)(z\, ■ ■ ■ , z n ) is 
holomorphic in D c . 

On the other hand, it is obvious that the recursive system of matrix Riccati equations and the 
system of independent matrix ODEs defined in proposition 2 depend analytically on Ai, • • • , A n G K 
and extend naturally to the case where Ai, • • • , A„, G C. Consequently, when (Ai, • • • , A n ) is replaced 
by {z\, • • ■ , z n ) G C n , the two systems determine solutions depending holomorphically on z±, ■ ■ ■ ,z n . 
It follows that tp(Xi, • • • , A n ) possesses a unique holomorphic extension to C™. 

Finally, since </>and tp coincide in the set {(Ai, • • • , A n ) G W 1 : Xj G (— c, c), j = 1, 2, • • • , n}, by 
the identity theorem, they coincide in D c . ■ 

With the preparations above, the proof of our main result on the formula for the joint char- 
acteristic function /(7 l5 • • • , 7 n ; Ai, • • • , A n ) defined in section one now follows easily. In fact, the 
result follows immediately from proposition 2 and lemma 3 with the observation that the set 

{OAi,-- - ,iA n ) : Aj = 1,2,- •• ,n} 

is contained in D c . 

3 An Example: the Two Dimensional Levy's Stochastic Area Pro- 
cess 

In this section, we are going to apply our result to study the two dimensional Levy's stochastic 
area process, first introduced by Levy in [10J. Namely, we consider the case where d = 2 and 

The Levy's stochastic area process is given by 

L t = [ {wPdWP - W^dWP), t ^ 0. 
J o 

We try to derive an explicit formula for the finite dimensional joint characteristic function of the 
coupled process {(Wt,L t ) : t ^ 0}. 

Throughout this section, = to < t\ < • • ■ < t n will be fixed. 

In section 2, the computation of the finite dimensional joint characteristic function of the coupled 
process {(Wt,Lf) : t ^ 0} reduces to the solution of a recursive system of symmetric Riccati 
equations and a system of independent first order linear matrix ODEs. In the case here, we will 
see that the Riccati system is actually real and scalar, and the linear system is explicitly solvable. 
In fact, we have: 
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Proposition 3.1 For the two dimensional Levy's stochastic area process {L t : t ^ 0} ; by using the 
same notation as in theorem 1 with the assumption 

A\ = • • • = A n = A, 

the solution matrices of the Riccati system are real diagonal matrices with identical diagonal entries, 
and the Riccati system essentially reduces to a system of real scalar Riccati equations recursively 
defined from j = n to j = 1 by 

d n n 

J t ki Aj ,..,iA n (t) = (A] + 2Aj A r)-2(Y, kiA r ,..,iA n (t))k iAj ,..,iA n (t) 

r=j+l r=j+l 

Moreover, the linear system in theorem 1 is explicitly solvable, namely, for j = 1, 2, • • • , n, 

ft n n 
H iAj ,..,iA n (t) = exp{ / (£K iAr ,... iiAn (s) + i(£ i A r )A)d8}, t € [0,tj]. 

Proof. We first consider the Riccati system. For j = n, 

CW(t)=( A ™ A °2 ), te[0,t n ], 
and the Riccati equation is defined by 

= ( A " ^)-a«-m(; - 1 ) 

-1 ) K ^(t))-Kf An (t), t€[0,t n ], 

K lAn (tn) = 0. 

By uniqueness, it is easy to see that the solution of the above equation is given by a matrix of the 
form 

where {h An (t) : t € [0,i n ]} solves the real scalar Riccati equation 

j t hA n {t) = A 2 n -kf An (t), te[0,t n ], 
ki An (t n ) = 0. 

For j = n — 1, by easy computation we have 

r t m _ ( A l + 2A n „iA n \ r , 

CiA^iAAt) - ^ Q A 2 + 2An _ iAn ), te [0,t n _i]. 
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It follows from uniqueness and the case j = n that -fQA n _i,iA„(i) is also a diagonal matrix of the 
form 

{t) = ( ^o A " ((> **..!«>)• te[o - 4 - 11 ' 

where {fciA„_i,iA n (*) : * G [0,i n _i]} solves the real scalar Riccati equation 

-^A n _ 1; iAn(*)' <G[0,*n-l], 
feiA„_i,iAn(*n-l) = 0. 

The rest of the argument follows from recursion easily (the crucial observation is that Ci\ jt ... ,iA n (t) 
is always a real diagonal matrix with identical diagonal entries). 

The second part of the proposition follows immediately from the fact that for j = 1, 2, • • • , n, if 
we denote as the coefficient matrix of the j— th linear ODE of the independent system, then 

m 

Now we are going to study the finite dimensional joint characteristic function 

fill,-- - ,7 n ; A ir-- ,A n ) 

of {(Wf, Lt) '■ t ^ 0}. To simplify our notation, we use {kj, Hj : j = 1, ■ ■ ■ , n} to denote the solution 
{kiAj,— ,«A n5 Hi\ jt ... ; jA„ : J = 1, • • • , in proposition 5. Assume first that (nondegeneracy) 

Aj + h A n 7^0. j = 1,... ,n. 

We first study the process {L t : i ^ 0}. A crucial observation is that for j = 1, • • • , n, by adding 
together from the j— th equation to the n— th equation in the scalar Riccati system in proposition 
5, we obtain a neat scalar Riccati equation without linear terms: 

Th Th Th 

!<!>(*)) = (I>) 2 - (Em*)) 2 , t € [o,^], 

r=j r=j r=j 

in which the unique solution is determined by the terminal data at t = tj. Let 

n n 

then it is not hard to derive that 

Cj sinh(cj(t — tj)) + Sj(tj) cosh(cj(i — tj)) 
3 3 Cj cosh(cj (t - tj)) + Sj (tj) sinh(cj (t — tj))' 

t e [0,tj],j = !,■■■ ,n, 
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where {sj(tj)}j =1 is defined recursively by 

8 u \ = c c i s inh (cj(*j-i ~ + s A t o) cosh(c j (t i _i - tj)) 
3 1 3 1 3 c j co^(c j (t j - 1 -t j ))+8 j (t j )w^(c j (t j - 1 -t j )Y 
j = 2,3, ■ ■ ■ ,n. 

Now apply theorem 1, we have 

n 

E[exp{^A fe L t J] 

k=l 

n 1 th 

= exp{^-/ Tr(K jAj ,.., lAn ( S ))^} 
j=l J o 

= ex P{y] / kj(s)ds} 

= expjy^ / Sj(u)du} 
j=i Jtj-i 

n 

= rr c i (*) 

fj[ cj cosh(c i (t j _i - tj)) + Sj(tj) sinh( Cj (i,_i - tj)) ' 

which seems, to our best knowledge, not appear in literatures, though it should follow from the 
Markov property and Levy's formula. 

It should be pointed out that the nondegeneracy assumption on {Aj : j = 1, ••• , n} is not 
important. In fact, if for some j, 

Cj = 0, 

and assume that {sj + i(t) : t £ [0, ij+i]} has been solved, then {sj(t) : t £ [0, tj]} can be sovled as 

Sj (t) = - 8 -fM - = - s -^) -, te[0,tj]. 

3 " 1 + Sj(tj)(t - tj) 1 + Sj +1 (tj)(t - tj) 

It is easy to see that {sj(t) : t £ [0, tj]} is actually the limit of the nondegenerate case as Cj — > 0. 
Moreover, the corresponding term in the product (*) can be written as 

1 



l + SjitjWj-i-tj)' 



which is also the limit of the nondegenerate case. Therefore, we still use the same notation even in 
degenerate cases. 

Now consider the coupled process {(Wt, L t ) : t ^ 0}. By proposition 5 and the above computa- 
tion, for j = 1, • • • , n, we can solve the linear system explicitly to obtain 
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where 



Oj(t) : 



/ Sj(u)di 
Jo 



ln c? cosh( Ci (t - tj)) + gjftj) sinh(cj(t - tj)) te < Qt . 
cjcoshicjtj) - Sj(tj)smh(cjtj) ' ' J 



Here the formula for {aj{t) : t £ [0, tj]} works in the degenerate case where Cj = as well. 

Finally, by theorem 1, we can write down the finite dimensional joint characteristic function of 
the coupled process {(Wt,L t ) : t ^ 0}. Namely, we have 

Theorem 3.2 For the coupled process {(Wt,L t ) : t ^ 0}, the finite dimensional joint characteristic 
function /(7 1; • • • , 7 n ; A 1 , ■ ■ ■ , A n ) is given by 

n 

fill,--' ,7n;Al,-- - ,An) = JJ 



?= i c i cosh(c i (t j _i - tj)) + Sj(tj) sinh( Cj (t^i - tj)) 
exp{-l f J {H*-\s)H*(t 3 )^ds}, 



where {cj, Sj(tj), Hj : j = 1, ■ ■ ■ ,n} is defined previously in this section, and {(ij : j = !,■■• , n} is 
defined recursively in terms of {7^, Hj : j = 1, • • • ,n} as in theorem 1. 
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